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ABSTRACT 

We investigate the effect of orientation-dependent selection effects on galaxy clustering 
in redshift space. It is found that if galaxies are aligned by large-scale tidal fields, then 
these selection effects give rise to a dependence of the observed galaxy density on the 
local tidal field, in addition to the well-known dependences on the matter density and 
radial velocity gradient. This alters the galaxy power spectrum in a way that is different 
for Fourier modes parallel to and perpendicular to the line of sight. These tidal galaxy 
alignments can thus mimic redshift space distortions, and thus result in a bias in the 
measurement of the velocity power spectrum. If galaxy orientations are affected only 
by the local tidal field, then the tidal alignment effect has exactly the same scale and 
angular dependence as the redshift space distortions in the linear regime, so it cannot 
be projected out or removed by masking small scales in the analysis. We consider 
several toy models of tidal alignments and orientation-dependent selection, normalize 
their free parameter (an amplitude) to recent observations, and find that they could 
bias the velocity amplitude f(z)G(z) by 5-10 per cent in some models, although most 
models give much smaller contamination. We conclude that tidal alignments may be a 
significant systematic error in redshift space distortion measurements that aim to test 
general relativity via the growth of large-scale structure. We briefly discuss possible 
mitigation strategies. 

Key words: large-scale structure of Universe - cosmology: theory. 



1 INTRODUCTION 



Redshift space distor tions (RSDs) have a long history as 
a cosmological probe. iKaiserl (|l987l ) first showed that even 
on very large scales, the galaxy power spectrum in red- 
shift space is significantly affected by bulk flows, and pro- 
posed this as a technique to measure fi m . The anisotropy o f 
the correlation function was explored by iHamiltonl i|l992l ). 
The redshift-space distorition technique was first applied 
to the Infrared Astronom i cal S atellite (IRAS) 2 Jy red- 
shift survey by IHamiltonl (| 19931 ). yielding a measurement 
of il m = 0.5^n'^K- Subs equent work with IRAS-selected sur- 



veys llCole et al. 19951) and optically-selected redshift sur- 
veys (|Lovedav et al.ll 19961 ; iRatcliffe et al.lll998l ) showed that 



the redshift-space anisotropy parameter f3 = Sl^ 6 jb was 
~ 0.5, providing strong evidence that either fi m < 1 or 
a high bias b ~ 2 was required. The 2-degree Field (2dF) 
survey represented a large st ep forward for redshi ft space 
distortion measurements, with lPeacock et"al] |200ll) report- 
ing j3 = 0.43 ± 0.07. Combining this approach with the 
galaxy bisp e ctrum , which independently gives b, allowed 
IVerde et all l|2002l ) to measure Q m = 0.27 ± 0.06, com- 
pletely internal to the 2dF survey. The Sloan Digital Sky 



Survey (SDSS) has reported RSD measurements from both 
magnitude-limi ted and color-selected (lumi nous red galaxy; 
LRG) samples (|Tegmark etalll2004 l2006j ) . RSD measure- 
ment s have also been reported for quasars (jda Angela et all 
l2005h and more recently (a nd with higher s i gnal-to-noise 
ratio ) for galaxies at z > 0.5 (Ross et al. 2 0071 ; IGuzzo et all 

boost ). 

Today, other methods of constraining f2 m are more 
powerful, namely combining the cosmic microwave back- 
ground (CMB) with super novae and baryon aco ustic oscilla- 
tion (BAO) measurements (jDunklev et al.|[2009l ). RSDs have 
also developed a reputation for having diflicult-to-control 
systematics since nonlinear evolution affects them at larger 
scales than the real space power spectrum. However, the dis- 
covery of dark energy has driven a revival in the subject of 
redshift space distortions. The same jumbo redshift surveys 
that are required by the BAO method will measure many 
modes of the galaxy field at large scales - precisely what 
is required to exploit the power of the redshift-space distor- 
tion techni que in the regime w here its systematics should be 



minimized (White et al 



2008). Modified gravity theories in 



which the acceleration of the universe reflects a breakdown 
of general relativity rather than a new contribution to the 
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Figure 1. The anisotropic selection effect. We show a Fourier mode k of the density field, oriented transverse to the line of sight in the 
left panel and radially in the right panel. The galaxies are shown aligned along the stretching direction of the tidal field, as appropriate 
for LRGs (but exaggerated). If selection effects prefer galaxies where the observer's line of sight looks down the long axis, then some 
galaxies whose short axes point toward the observer are lost (marked with an X in the figure). The remaining galaxies show stronger 
clustering for transverse modes (fc _L ft) than for radial modes (fc || n), thereby contamination redshift-spacc distortion measurements. 
In this case the effect is opposite to the linear Kaiser effect and hence biases estimates of the velocity power spectrum low. 



cosmic energy budget have motivated studies of the growth 
of structure, as it is impossible to distinguish these mod- 
els from da r k ene r gy using distance information alone (e.g. 
Ishak et al . 2006; Song fc Percival 2008; Song fe Kovamal 



20091 ). In this case redshift-space distortions generally mea- 



sure the velocity power spectrum P v (k). Finally, advances in 
understanding nonlinear evolution with lar ger iV-body sim- 
ulatio ns and with halo occupation models ( Coorav fc Shethl 
2002), and novel techniques using multiple galaxy tracers to 
reduce the uncertainty on P v (k) (|McDonald fc Seliakll2008l ) 
have made the case for RSDs more compelling. Indeed, RSDs 
are part of the prime science case of the proposed Euclid 
dark energy mission. 

RSDs suffer from a very different set of limitations than 
alternative probes of the growth of cosmic structure such as 
cosmic shear (CS), the integrated Sachs- Wolfe (ISW) effect, 
and the cluster mass function. CS directly probes the matter 
distribution without assumptions about galaxy biasing; it is 
however technically difficult to measure galaxy shapes with 
sufficient control of systematics, and practical implementa- 
tions of CS rely on photometric re dshifts. The ISW effect has 
turned out to be relatively clean (|Giannantonio et al,ll2008l ; 
IHo et al"1|20 08i), but its statistical power suffers from severe 
cosmic variance limitations. Clusters are easier to observe 
than CS, and are numerous so that statistical errors can be 
made small even with present data; the main uncertainty is 
instead in the astrophysics of relating cluster observables to 
virial masses. 

The most worrying systematic error in the RSD method 
is the behaviour of galaxy biasing and velocities in the non- 
linear regime. However, we show in this paper that the align- 
ment of galaxies by large-scale tidal fields results in an addi- 
tional systematic error. If galaxies are preferentially aligned 
along the stretching axis of the tidal field, as is the case 
for LRGs, and there is any viewing direction-dependent se- 



lection effect (e.g. preferring galaxies where we look down 
the long axis, as would occur with an aperture magnitude 
cut), then Fourier modes of the density field along the line 
of sight are suppressed because the galaxies in the troughs 
are more likely to be selected than galaxies in the crests. 
This is shown in cartoon form in Fig. ff] Tidal alignments 
of disk galaxies would if present produce a similar effect be- 
cause these galaxies are often selected by optical continuum 
or emission line luminosity cuts, which suffer inclination- 
dependent internal extinction. The alignment effect is small 
- we will argue that in pessimistic scenarios it contaminates 
redshift space distortion measurements by 5-10 per cent - 
but we will show that it exactly mimics the angular and scale 
dependence of RSDs, making it hard to remove. Moreover, it 
is not included in mock catalogues that paint galaxies onto 
JV-body simulations. 

This paper is organized as follows. In Section [2] we 
outline the standard theory of the RSD technique. In Sec- 
tion [3] we show how they can be altered by galaxies with 
tidal alignments and viewing direction-dependent selection 
effects. Section|4]specializes to the physically motivated case 
of galaxies aligned by the large-scale tidal field. In Section 
we consider some crude models of orientation-dependent se- 
lection effects and find that they may be significant. Sec- 
tion [5] compares the effects of tidal alignments on RSD to 
their more familiar effect on CS. We conclude and discuss 
mitigation strategies in Section [7] 

The early parts of this paper, especially Section [3] are 
intended to be general, exploring all models allowed by sym- 
metry; they are unavoidably heavy on mathematical for- 
malism. Section [5] in contrast, is intended to consider spe- 
cific examples; it uses rough calculations since orientation- 
dependent selection effects are not amenable to precise cal- 
culation. Section [6] is included mainly to develop intuition 
and hence also makes use of some rough calculations. 
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2 STANDARD REDSHIFT SPACE 
DISTORTIONS 

We begin by revi ewing the st andard theory of redshift- 
space distortions l|Kaiserl Il987t ). In this paper, we denote 
the galaxy bias by b and the rate of growth of structure by 
/ = dlnG/dlna, where G is the growth function. 

We consider an underlying matter density fluctuation 
field 5 m (x), where x — (xi,X2,X3) is the comoving posi- 
tion. We work in the flat-sky approximation and take the 
observer's line of sight to be in the a^-direction. The bulk 
velocity field is described in linear perturbation theory by 

v(x) = -/V.V« 2 J m (x). (1) 

For a sample of galaxies with linear bias b, the real- 
space density fluctuation of galaxies in the linear regime is 
5g (as) = bdmlx). However, in redshift space, the density is 
also corrected by the Jacobian of the conversion from real 
to redshift space, 1 — <?i;3 /<9:r3 Q That is, 

* g (x) = b5 m (x) - (2) 
0x3 

These equations are most easily represented in Fourier space: 
v(k)=-if±6 m (x) (3) 

and 

5 g (fe) = M m (fe) - ik 3 v 3 {k). (4) 

The overall galaxy density in redshift space is then 

i g (fe) = (b + / M 2 )Mfc), (5) 

where \i = ks/k is the cosine of the angle between the line 
of sight and the direction of the Fourier mode under consid- 
eration. This implies an observed galaxy power spectrum, 

P g (fe) = (6 + / M 2 ) 2 P m (fc). (6) 

The amount of anisotropy depends on the parameter /3 = 
f/b. 

The standard use of redshift space distortions is to use 
the amplitudes A 2 (fc) = fc 3 P g (fc)/27r 2 , which satisfy 

Ag(fe) = 6A m (fc) + m 2 / A m (fc), (7) 

to extract the quantity /A m (fc). If the matter power spec- 
trum is normalized at high redshift by the CMB, then 
A m (fc) oc G, and redshift space distortions can be used to 
obtain the quantity f(z)G(z). Any contamination to Eq. (J7J) 
produces a fractional error in f(z)G(z) equal to the frac- 
tional change in the /x 2 term, to the change in the coefficient 
of fi 2 divided by /A m (fc). 

The above analysis assumes that when we measure 
a Fourier mode of the galaxy density field we know the 
calibration of the radial and transverse distance scales so 
that the comoving k is known. This could be done using 
the BAO combined with a ruler length set by the CMB, 
or in a parameterized cosmology by combining all cosmo- 
logical probes. The parameter f3 only requires the relative 
calibration of k« and k±, which amounts to determining 
the product H(z)D(z) where H(z) is the Hubble rate and 

1 In linear theory we may neglect the change in the position x 
as we move from real to redshift space. 



D(z) is the comoving angular diameter distance. In princi- 
ple this could be done int ernally to the galaxies using the 
lAlcock fc Paczvnskil (|l979l ) test. In this case, extracting / 
would require independent information on the galaxy bias, 
e.g. from the three-point function. 



3 THE PROBLEM 
3.1 Intrinsic alignments 

If galaxies were randomly oriented, Eq. © would fully 
describe the linear regime. However, some classes of real 
galaxies, including LRGs, exhibit correl ations of their ori- 
entation with larg e-scale structure (e.g. lBinggelilll982l ; see 
iHirata et al.l 120071 for a recent measurement of the correla- 
tion function). This section defines the notation needed to 
handle selection effects that depend on galaxy orientation 
and viewing direction. 

The galaxy orientation is described by 3 Euler angles 
(6,(j),ip). These define a rotation matrix Q(0,cj>,ip) £SO(3) 
that transforms the "lab" frame coordinates to a frame 
aligned with the galaxy; see Appendix [D] for explicit ex- 
pressions. The well-known volume element of SO (3) is 

d 3 Q = sin^d^d^d^, (8) 

and its total volume is f so , 3 -, d 3 Q = 87r 2 . We consider the 
conditional probability distribution V(Q\x) for the orienta- 
tion Q of a galaxy at position x. In the case of randomly 
oriented galaxies we have P(Q\x) = 1/(8tt 2 ), but m gen- 
eral this probability distribution may depend on the local 
environment (tidal field, density, etc.) at x. 

Our next ingredient is the dependence of the observa- 
tional selection function on the viewing geometry and galaxy 
orientation. We assume that the observer is looking along 
the line of sigh10 n. The appearance of the galaxy then de- 
pends on the observer's line-of-sight expressed in the galaxy 
frame, i.e. m = Qn. It is therefore possible for the proba- 
bility of selecting a galaxy to vary depending on Qn: 

P <xl + T{Qh,x). (9) 

We may take the anisotropy of the selection function to have 
mean zero when averaged over all possible viewing direc- 
tions, 

/ T(m,x)d 2 m = 0, (10) 

since any m-independent constant added to T can be ab- 
sorbed into the local galaxy density. Note that Eq. (|10p im- 
plies that the selection probability variation averaged over 
all possible orientations is also zero, i.e. 

[ T(Qn,a;)d 3 Q = 0. (11) 

JsO(3) 

The tidal alignment effect on large scale structure ob- 
servations depends on the combination of the intrinsic align- 
ment model, which determines V(Q\x), and the dependence 
T of the selection function on geometry. The number density 

2 We define this so that the observer points their telescope in 
direction n; the light from the galaxy is propagating in direction 
— n. 
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of selected galaxies TV(selected) is then related to the true 
number density of galaxies JV(true) by averaging Eq. ((9} 
over the distribution of galaxy orientations at x: 

N (selected) 



7v(true) ^ / sow nQWIl + r(Q A>1B )]d»Q 

= 1+ / P(Qjx)T(Qn,:r)d 3 Q. (12) 

JSO(3) 



This motivates defining 



e(n\x) = / V(Q\x)T(Qn,x)d- , Q. 

JS0(3) 



(13) 



which represents the viewing direction dependent selection 
function at position x. Examining Eq. (lll[l . we see that 
e vanishes if either the galaxy orientations are isotropi- 
cally distributed [in which case P(Q|x) = 1/(8tt 2 )], or if 
T vanishes. Both intrinsic alignments (V depends on Q) 
and orientation-dependent selection (T ^ 0) are required to 
produce an effect on large scale structure. 

[In the case of an axisymmetric disk galaxy, the third 
Euler angle ip is irrelevant. In this case, integrals such as 
Eq. (|13fl reduce to integrals over the disk normal vector L.] 

The real-space density of observed galaxies at position x 
as seen by an observer looking in direction n has fluctuation 



l + S^ ohB \x,n) = [l + ««(as)][l + e(ft|aO]. 



(14) 



In a real galaxy survey, the actual viewing direction n is ez, 
and our final computation of the galaxy properties will be 
obtained by substituting n = e$ into Eq. (|14p. 

Equation (I10|l implies that e(n\x) to average to zero 
over the unit sphere (i.e. over all possible viewing directions) , 



e(n\x)d n = 0. 



(15) 



We then find, instead of Eq. ([5]), the observed redshift- 
space galaxy density in the linear regime is 



tf g (fc) = (6 + /^ 2 )<5 m (fe) + e(e 3 |fe) 



3.2 Statistics of intrinsic alignments 



(16) 



In order to construct the observed galaxy power spectrum, 
we need to develop the statistics of e(n\x). This is a random 
field as a function of a position x and direction n, and its sta- 
tistical description can be achieved in the same way that one 
decomposes the CM B temperature 8T(Xj n) in cosmological 
pertu rbation theory (|Ma fc Bertschinger|[l995l ; Irlu fc White! 
1 19971 ). One first Fourier transforms the spatial variables to 
get e(n\k). Then one introduces a rotation matrix R GSO(3) 
that rotates the unit vector k to the e-j, axis, i.e. Rfc = e.j,\2 
Here R depends on k, but we will not write this dependence 
explicitly to avoid confusion. Finally, we decompose the de- 
pendence of e(njfe) on the viewing direction n into spherical 
harmonics, 



e~(n|fc)=£ £(- 



Z — 1 m — — l 



4tt 
21 + 1 



e tm (k)Yi m (Rn). (17) 



3 Note that R is not unique, since it amounts to an arbitrary 
choice of which unit vector in the plane orthogonal to k gets 
rotated to e\. 



This corresponds to choosing the "North Pole" of the spher- 
ical harmonic basis to be in the direction k. We note that 
Eq. (|15|l eliminates the I — term. 

Translation invariance forces different fc-modes to be 
independent, and rotational invariance around k forces dif- 
ferent different values of m to be independent. Therefore, 
we can write the power spectrum of e: 

(er m (fc)e !W (fc')> = (27r) 3 P e iiV (fc)<5 mm ,5 (3) (fc - k'). (18) 

These power spectra are not arbitrary. Appendix |X] shows 
that since e(n|x) is a real field, P" m (k) must also be real. 
Invariance under reflection across a plane containing k forces 
the restriction 



P l / m {k) = P l /'' m (k), 



(19) 

and by swapping the I and labels in Eq. (jTHJ and recall- 
ing that the power spectra are real, we find the symmetry 
relation: 



Pf' m (k) = P l c ' lm 



(k). 



(20) 



Therefore the power spectra are completely described by the 
cases with ^ m ^ I' ^ I, with I, I' > 0. 

In addition to its intrinsic power spectrum, e can have 
a cross-power spectrum with the matter field. This can be 
defined in analogy to Eq. (jTHJ ■ The matter field has the same 
symmetry as the I = m = mode of e (if it existed), so by 
rotational symmetry only the m = components of e can 
correlate with the matter field. We thus have 

(&(fc)«o(fc')> = (2^) 3 P,L(fc)5 (3) (fe - *') (21) 

for I > 0. Again, P^(k) must be real. 

A further simplification is possible for galaxies that are 
inversion-symmetric, i.e. have the same appearance if viewed 
from the opposite direction. This is true for optically thin, 
triaxial elliptical galaxies, and is true for disk galaxy models 
that are axisymmetric and have a reflection symmetry across 
their equatorial plane. It may also be true in a statistical 
sense for real disk galaxies: while the specific arrangement 
of dust and H II regions will not have inversion symmetry, it 
is plausible that these local imperfections would be unique 
to each galaxy and not correlated with large-scale structure. 
In this case, e(ri|x) = e(—n\x) and hence e; m (fe) = for 
odd I. All power spectra with odd I or I' then vanish. 

In traditional cosmology language, the m — contri- 
butions to this equation are scalars, m = ±1 are vectors, 
m = ±2 are tensors, and so on. One would expect in the 
linear regime that only the scalars (m = 0) are present 
since large-scale density perturbations are scalars and can- 
not source any alignment of galaxies with m/0 symmetry. 
(Physically, in the absence of mode coupling, a Fourier mode 
k possesses rotational symmetry around the axis k that pre- 
vents any m ^ component.) At small scales it is possible 
that the m 7^ alignments arise from nonlinear mode cou- 
pling. Since Eq. (|16|) is valid only when the density and ve- 
locity fields are in the linear regime, in the body of this paper 
we will only consider the scalar (m = 0) intrinsic alignments. 
We note, however, that the scale at which nonlinear effects 
start to influence intrinsic alignments is unknown and may 
be larger than the traditional t m !„ ~ 10/i -1 Mpc used in 
galaxy surveys. In this case one might want the equations de- 
scribing a general intrinsic alignment model; these are given 
in Appendix [B] 
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3.3 Galaxy power spectrum 

Having defined the intrinsic alignment power spectra, we 
now compute the observed redshift-space power spectrum 
of the galaxies, i.e. of Eq. (|16[1 . We consider here only scalar 
intrinsic alignments with even / only (i.e. we assume inver- 
sion symmetry). We first plug in Eq. (I17|l . 



Sg(fc) = (6 + /M 2 )5 m (fe) + J2 ^o(fe)Fw(Re 3 



(22) 



1^2 2, even 

Taking the power spectrum gives 
P g (fc) = (6 + / M 2 ) 2 P m (fc) 

+2(6 + / M 2 )5R ]T (- 



i) i J^^Pi E (fc)F !0 (Re 3 ) 



2Z + 1 



47rP e '''°(fc)y,g(Re 3 )y,, (Re 3 ) 
a/(21 + 1)(2Z' + 1) 



(23) 



Here the first term is the conventional redshift-space power 
spectrum, the second comes from correlations of the matter 
(and velocity) fields with intrinsic alignments, and the final 
term is the pure intrinsic alignment term. 

Equation (I23|l can be simplified by replacing Y;o with a 
Legendre polynomial^ We can evaluate its argument via 



IX = e 3 ■ k = (Re 3 ) ■ (Rfc) = (Re 3 ) • e 3 , 



(24) 



where the first equality is the definition of fi, the second 
is the definition of a rotation matrix, and the third follows 
from the condition Rfc = e 3 . Thus 



4tt 



r i0 (Re 3 ) = Pi(n). 



21 + 1 

This substitution gives rise to 
P g (fc) = (6 + / M 2 ) 2 P m (fc) 

+2(& + / M 2 ) J2 (-i) l/2 PLMPi(ri 

1^2, even 



(25) 



J2 (-if-^P^^Pti^itx). (26) 



4 TIDAL ALIGNMENT MODELS 

The best-motivated model for intrinsic alignments is to sup- 
pose that large-scale tidal fields would induce a preferential 
direction in galaxy formation. The selection probability for 
a galaxy would then depend on the both n and the con- 
figuration of the tidal field. Tidal fields from linear regime 
density fluctuations are small: by definition to be linear the 
tidal field is ^ H 2 where H is the Hubble constant, and 
hence must also be less than £~T 2 where t<j is the relevant dy- 
namical timescale of the galaxy at any stage of its collapse. 
Moreover, the small-scale nonlinear tidal fields are stronger 
than the large-scale tidal fields because fc 3 P m (fc) is a rapidly 
increasing function of A:. A plausible model for tidal align- 
ments would then be to take the large-scale tidal field as 



a perturbation and Taylor expand to lowest nonvanishing 
order: 



e(h\x) 



A 



AirGa 2 p m (a) 



hihjViV j 



9(x), (27) 



where 3 , (x) is the Newtonian gravitational potential, A is 
the expansion coefficient, and the denominator 47rGa 2 p m 
serves to make A dimensionless0 The second derivative of 
is the tidal field tensor, and the specified contraction with 
h is the only possibility consistent with the vanishing angle- 
average of e(n|x) since the tidal field is a quadrupole. Note 
that even for the same population of galaxies, it is possible 
for A to be redshift-dependent. This dependence cannot be 
predicted without a detailed microphysical model of galaxy 
alignments. Using the Poisson equation, Eq. (|27[1 can also 
be re-expressed as 



e(n\x) = Asij(x)riinj 



where 



s l3 = ( V^V- - -S i3 ) 5 m (x) 



(28) 
(29) 



is the dimensionless tidal field. 

An intrinsic alignment model of the form of Eq. (|27[1 
could also be motivated b y renormalization group arguments 
(|McDonald fc R ov 2009): at first order in perturbation the- 
ory, the density 8 and tidal field tensor Sij fully describe the 
history of any patch of material in the Universe. Presum- 
ably second and third-order terms could be incorporated 
into Eq. (|27[1 and used to define a renormalized intrinsic 
alignment model in analogy to the [McDonal d fc Rovl (120091 ) 
approach to the galaxy density. We note that for disk galax- 
ies, an inclination-dependent selectio n effect combined with 
the quadratic tidal a l ignment model (tCrittenden et al.ll200ll : 
ICatelan et"afl 1200 ll : iMackev et ail l2002h is an example of 
such a second-order term. Like the second-order biasing of 
galaxies, quadratic alignment is most prominent on small 
(quasilinear) scales. We briefly return to this issue in Sec- 
tion [5]2] A full analysis of the quasilinear regime, including 
quadratic alignments and their effects on e.g. the galaxy bis- 
pectrum and quasilinear RSDs, is beyond the scope of this 
paper. 

We now estimate the effect of tidal alignments on the 
observed galaxy power spectrum. In Fourier space, Eq. (|28|l 
becomes 



e(n|fe) = A 



(A-k) 2 -^ 



5 m (k) = -AP 2 (n-k)S m (k). 



(30) 



Expanding in terms of spherical harmonics, this implies 



e 20 (fe) = --A(5 m (fe), 



(31) 



with all other components equal to zero. The implied power 
spectra are 



PL(k) = ~AP m (k)S l2 



(32) 



4 Following convention, we denote the Legendre polynomials by 
Pi and associated Legendre polynomials by P ; m . By context these 
should not be confused with power spectra. 



5 Whether in practice a Taylor expansions such as this is appro- 
priate for all types of galaxies is an open question. For example, 
higher-order terms in the tidal field could dominate, or the actual 
tidal field could lie outside the radius of convergence of the series. 
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and 



Pf' m (k) = ^A 2 P m {k)8 l2 8 ll2 8 m0 . 



(33) 



The fact that only / = 2 components appear is a direct result 
of the alignment model being linear in the tidal field, which 
transforms as a quadrupole. Substitution into Eq. (|26[1 gives 



P g (fe) = (6 + //i 2 ) 2 P m (fc) + ^(6 + / M 2 )P 2 ( M )P m (fc) 



+lA 2 [P 2 ^)fP m {k). 



(34) 



Because e is completely determined by <5 m in this model, this 
simplifies to 



b-^ + (f + A)» 2 



Pm(fe). 



(35) 



We conclude that in the tidal alignment model, the ob- 
served redshift-space power spectrum retains the same func- 
tional form, but the intrinsic alignments alter the coeffi- 
cients. The "obvious" systematics tests for redshift space 
distortion measurements, such as testing the scale depen- 
dence of f3 or the "constant+p 2 " dependence of sj P g (k), will 
not detect this effect. The intrinsic alignments change the 
normalization of the real-space (p = 0) galaxy power spec- 
trum, but this is degenerate with bias. They also change 
the apparent rate of growth of structure, /[P m (fc)] 1//2 — > 
(/ + A)[P m (k)}^ 2 , which is a real systematic error whose 
fractional magnitude is |^4|//. One can also note the sign 
of the effect: if galaxies are preferentially selected in regions 
where the tidal field is compressional along the line of sight 
(A > 0) then the RSD is enhanced; whereas if galaxies 
are preferentially selected in regions where the tidal field 
is stretching along the line of sight (A < 0) then the RSD is 
suppressed. 



5 MODELS FOR A 

Now that we have established the effect of tidal alignments 
on the galaxy power spectrum, it is time to consider the 
likely value of the A parameter. Since this parameter de- 
scribes how tidal fields feed into galaxy orientations and 
ultimately observational selection criteria, we consider two 
separate cases. The first is an "elliptical galaxy" case in 
which selection criteria are dependent on e.g. isophotal mag- 
nitudes or effective radius, which appear different depending 
on whether one looks down the long axis or the short axis of 
the elliptical. The second is a "disk galaxy" case in which the 
selection criteria are blue continuum or emission line inten- 
sities that suffer inclination-dependent extinction. In both 
models, we can roughly predict the normalization of how 
much galaxy orientations affect selection but not how much 
tidal fields affect orientation; for the latter, we take observa- 
tional results on the correlations of galaxy ellipticities with 
large scale structure. 

Note that all of the models considered here are merely 
examples and may not represent the actual contamination in 
a particular survey. The results depend sufficiently strongly 
on the detailed selection criteria that each survey must do 



its own calculation, either along the lines of those presented 
here or using simulations. 



5.1 Elliptical galaxy 

We model an elliptical galaxy as an optically thin (i.e. neg- 
ligible dust) triaxial system. The volume emissivity at a po- 
sition s relative to the centre of the galaxy is taken to be 
j(s) = J~{p), where the ellipsoidal radius p is p 2 = s ■ V _1 s. 
The matrix V is symmetric and unimodular (det V = I) and 
contains information about the anisotropy of the galaxy, 
whereas J(p) specifies the radial profile. For a spherical 
galaxy, V = I. We define the deviation from spherical sym- 
metry to be W = V I , and will work in the limit where W 
is small. To first order, W is traceless. It is described by 5 
numbers: 2 linearly independent eigenvalues that specify the 
amount and type (oblate vs. prolate) of triaxiality, and 3 Eu- 
ler angles. We will argue here that for small W, the selection 
probability T(Qe3, x) is proportional to W33(x). Then in ac- 
cordance with Eq. (|13[) we will find that e(e3|x) oc (W33)(x), 
where the coefficient of proportionality depends on the se- 
lection algorithm. We will argue that (Wij)(x) oc Sij(x) (see 
Appendix [D] for an explicit demonstration in terms of the 
integral over Q). The coefficient of proportionality, needed 
to determine A, will be fixed using measurements of galaxy 
ellipticity-density correlations. 

We do not observe the full 3-dimensional structure of 
the galaxy, but rather its 2-dimensional projection 



Hs±) 



j{s± + s 3 e 3 ) ds 



3, 



(36) 



where s± is a position in the plane of the sky, i.e. the 12- 
plane. 

In Appendix [C] we show that the observed two- 
dimensional projection is 



I{a±) = \1+ -W33) /o(c). (37) 

Here 

7 (c) = J J(\c + s 3 e 3 \) ds 3 (38) 
is the projected image of the spherical galaxy model, 



i - §W ± ) si (39) 



is a rescaled skewer position in the 12-plane, and Wj_ is the 
projection of W into the 12-plane (i.e. the 2x2 submatrix). 
Note that Wi need not be traceless; indeed TrWi = — W33. 

The observed properties of this galaxy are as follows. 
The effective radius r c is reduced from its spherical version 
in accordance with Eq. (|39p: 



/■ = 1 1 - ^ it'i.i j /■,.!,. 



(40) 



The total flux F of the galaxy is unchanged. Finally, the 
ellipticities of the isophotes (surfaces of constant c) compo- 
nents are 

ei = ~(Wn - W22) and e 2 = Wyt. (41) 

LRGs are usually selected by a combination of cuts in 
colour space, combined with a magnitude cut. As long as 
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the galaxy is optically thin, there is likely to be very little 
orientation-dependent change in colour. However, the ap- 
parent magnitude of a galaxy may depend strongly on its 
orientation. Here we provide some examples. 

Model magnitudes: If galaxies are selected based on 
model magnitudes, and the model fit accurately describes 
the radial profile, then the model-fit magnitude corresponds 
to the orientation-independent flux F. This gives e = 0. 

Petrosian magnitudes: For galaxies described by a sin- 
gle radial profile J(p) (e.g. a deprojected de Vaucoulcurs 
model), that are well-resolved (in the sense that the Pet- 
rosian radius is well beyond the smearing of the central cusp 
due to the point-spread function), and have small ellipticity 
(W small), the Petrosian magnitude should be invariant un- 
der the rescaling of r e described by Eq. (|40|l . This is because 
the Petrosian radius will be a fixed multiple of r e and will 
capture a fixed fraction of the galaxy's flux. Within these 
approximations wc would have e = 0, but the robustness 
of this conclusion should be tested in simulations if high 
accuracy is required. 

Aperture magnitudes: If galaxies are selected based on 
aperture magnitudes, then we count more of a galaxy's light 
when it is viewed down a long axis than a short axis. As 
an example, suppose that the aperture is cut off at y effec- 
tive radii the spherical galaxy model, and the shape of the 
galaxy's profile is such that the fraction of the luminosity 
within yr c is Q(y)- A 1 per cent decrease in r e results in a 
1 per cent increase in y and hence a yQ' (y)/ Q(y) per cent 
increase in Q(y). Then using Eq. (|40|l . the aperture flux F ap 
is modified according to 



apO 



4 Q(y) 



(42) 



If the galaxies are selected according to a flux limit and the 
cumulative luminosity function of galaxies has logarithmic 
slope rj, i.e. dlnn/dlnF m i n = —t], then a 1 per cent increase 
in measured flux for galaxies near threshold corresponds to 
an rj per cent increase in the number density. Thus 



4 Q(y) {W33> 



(43) 



Isophotal magnitudes: Isophotal magnitudes are similar 
to aperture magnitudes in that we count more of a galaxy's 
light when it is viewed down a long axis than a short axis. 
However, in this case the cutoff for flux determination de- 
pends on the surface brightness profile. The surface bright- 
ness at y effective radii is 



F Q'(y) 



V 



(44) 



If we change r e but hold the surface brightness constant (as 
appropriate for determining how the isophotes move), we 
find 







SI (y) = _ 2 — + 
I(y) r c 



dy y 



Sy. 



(45) 



The quantity in brackets evaluates to Q" (y) / Q' (y) — y~ ■ 
We also know from Eq. (|40|) that Sr c /r c = — W33/4. Thus 
we may solve for 5y: 



5y 



-W, 



2[Q"(y)/Q'(y) - y~ 



(46) 



Selection factors for de Vaucouleurs profile 




Figure 2. The selection factors \ of Eq. 1491 1 for a de Vaucouleurs 
galaxy. Either aperture or isophotal magnitudes are used for the 
flux cut. The factor depends on y, the number of effective radii at 
which the photometry is cut off for a typical galaxy at the faint 
end of the sample. If model or Petrosian magnitudes were used, 
we would have \ = 0. 



This implies a fractional change in isophotal flux of 
F is Q'(y)Sy Q'(y)W 33 

F is0 Q(y) 2Q{y)[Q"{y)/Q'(y)-y- 1 ] 

For a luminosity function slope rj, we find 
riQ'iv) 



(47) 



—AW 33 ). (48) 



2Q(y)[Q"(y)/Q'(y)-y 
Note that the coefficient of W33 is positive, because 
Q" (y) / Q' (y) ~ V 1 < for any surface brightness profile 
that decreases as one moves outward (e.g. de Vaucouleurs). 
In all of these cases, one may write 



VXiWas), 



(49) 



where \ depends on the method of computing fluxes (model, 
Petrosian, aperture, or isophotal), the galaxy profile, and 
(in the latter two cases) the typical number of effective 
radii y at which the flux is computed for galaxies near 
the threshold. For the lde Vaucouleursl (Il948l ) profile, where 
Q'(y) oc y exp(— 7.67j/ 1//4 ), values of \ are shown in Fig- 
ure [2] Note, however, that Petrosian, aperture, or isophotal 
magnitudes of galaxies at radii affected by the point-spread 
function of the telescope would exhibit more complicated 
behavior that depends also on the dimensionless resolution 
factor (model magnitudes will be affected too if the model 
does not correctly describe the galaxy). 

We now suppose that the mean value of W for galax- 
ies in a particular region of space has some dependence on 
the tidal field surrounding it. Since both Wij and the tidal 
field Sij are traceless-symmetric tensors and have quadrupo- 
lar symmetry, the lowest-order allowed term in the Taylor 
expansion of Wij as a function of Sij is: 



(Wij) = 2Bsij = 2B ( ViVjV"^ - - ) 5 m (x). (50) 



An alternative argument for Eq. (|50[) in terms of the Euler 
angles can be found in Appendix [D] We also argue there 
that although the full specification of the alignment of a 
triaxial galaxy to linear order in Sij requires 2 parameters, 
the statistical effect on (Wij) only involves 1 parameter (B). 
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Comparing Eq. flSD) to Eqs. (J28J) and (g9j) gives the 
proportionality coefficient, 

A = 2 VX B. (51) 

As discussed in Sec. 13.11 there is an effect on large scale 
structure only when the galaxies have nonrandom orienta- 
tions (B 7^ 0) and when there is an orientation-dependent 
selection effect (x ^ 0). 

The normalization constant B depends on the details 
of elliptical galaxy formation. Rather than choosing a value 
based on theoretical considerations, we normalize B to mea- 
surements of intrinsic ellipticity correlations of elliptical 
galaxies. Using Eq. (|41[) . we find that for a Fourier mode 
k perpendicular to the line of sight, e.g. in the 1-direction 
(ft = feei), we have ei(fc) = i?<5 m (fc). This implies a cross- 
power spectrum between matter and ellipticity of 

P me (fe) = BP m (k). (52) 

This power spectrum - or more prec i sely, the cross- 
correlation - was measured bv lHirata et al.l ((2007) , who were 
interested in intrinsic alignments as a contam inant of cos- 
mic sh ear measurements. For the latter reason, iHirata et al.l 
l|2007h presented not the correlation with ellipticity compo- 
nents ei,2 but with "intrinsic shear" 7 lj2 = ei, 2/1-74, where 
1.74 is the shear responsivity. Note that in the notation of 
iBernsteinl (|200ST ). b K = B/1.7A. The projected correlation 
function ws+(r p ) of the matter and intrinsic shear is 

w s+ (r p ) = ~ J P m (k)J 2 (kr p )kdk. (53) 

IHirata et al.l l)2007h showed that the scale dependence im- 
plied by Eq. (|53|l is a good fit to the intrinsic alignments 
of LRGs. They normalized their correlation function by its 
value at r p = 20h~ 1 Mpc and z = 0.3, and usin g the non- 
linear matter power spectrum (Sm ith et al 1 l2003h . Eq. ((531) 
gives a correlation function of — 3.3b K /i _1 Mpc. Comparing 
to the observed value of 0.059(L/L ) 1 ' 48 ^" 1 Mpc (where L 
is the luminosity K + e-corrected to z = and Lo corre- 
sponds to a corrected r-band absolute mag nitude M°° of 
— 22) suggests that 

/ r \ 1.48±0.64 

b K = (-0.018 ± 0.006) (— J . (54) 

Thus at z = 0.3, B w -0.03 for LRGs of magnitude M° ° = 
-22, rising to -0.06 at M°'° = -22.5 and -0.12 at M° ° = 
-23. 

In order to estimate the anisotropic selection param- 
eter A, we additionally need the slope of the cumulative 
luminosity function r\ — — dlnw /dln Fmin- We ma y use 
the LRG luminosity function of IWake et al.l (|2006T lFl At 
Mj'° = —22.5, we find 77 = 4.0. If galaxies are selected 
with isophotal magnitudes with a surface brightness cutoff 
at ~ 3r c , then \ ~ 0.05. Combining with B ~ —0.06, this 

leads to A' 0.024. Since at z = 0.3, / « 0.65, we find that 

the fractional contamination to the redshift-space distortion 

6 IWake et al] 11200611 note that for typical L RGs, M° 2 M° ° + 
0.11; but note that the absolute m agnitude in IHirata et ail (1200711 
is referenced to lO/i" 1 pc, whereas lWake et al.l d2006T) refers abso- 
lute magnitudes to 10 pc assuming h = 0.7. This intro duces an ad- 
dition al offset of 5 log 10 0.7 so that the magnitudes in lWake et al.l 
j2006h are given by M^ ake = M®'° — 0.66. 



measurement \A\\/f is ~ 4 per cent. For LRGs a half a mag- 
nitude brighter at M°'° = -23, 77 = 4.5 and B ' 0.12, 

so the contamination increases to ~ 8 per cent. If one used 
aperture magnitudes instead, with a typical aperture cutoff 
at ~ 3r c , then we would have x ~ 0-06 and our estimated 
fractional contamination \A\/f for these two cases would rise 
to 5 and 10 per cent, respectively. 

These values are extremely rough and are merely in- 
tended to show that there exist regimes in which the cosmo- 
logical interpretation of the redshift-space distortions is sig- 
nificantly altered. The level of contamination for a particular- 
survey depends in great detail on its selection criteria; if it is 
significant, r/x would best be determined by a combination of 
Monte Carlo simulations and direct use of the observed dis- 
tribution of galaxies in colour-magnitude space rather than 
the simplified analytic arguments used here. Also note the 
sign: for isophotal or aperture magnitude selection, where 
galaxies are more likely to pass cuts if aligned in the radial 
direction, we have A < 0. 

The leading spectr oscopic LRG surveys t oday, such as 
the SPSS LRG surv ey l|Eisenstein et al] 1200 il l and 2SLAQ 
(|Cannon et al] [2006). use Petrosian or model magnitudes 
for their selection. This dramatically reduces the intrinsic 
alignment contamination since the selection should be al- 
most orientation-independent (except for discrepancies be- 
tween the model and the actual galaxy profile, and for the 
handful of objects with significant dust). 

5.2 Disk galaxy 

We now consider an alternative model for anisotropic selec- 
tion: a disk galaxy selected in either emission lines (e.g. [O 
n] or Ha) or optical/ultraviolet continuum will suffer less 
extinction if viewed face-on. Therefore selection based on 
the apparent magnitude or emission line flux favors face-on 
galaxies. The flux from the galaxy depends on the detailed 
distribution of stars or star-forming regions and dust. For 
some applications the extinction and scattering problems 
can be reduced by working in the rest-frame infrared, how- 
ever in large redshift surveys that require blue or ultraviolet 
conti nuum as part of the colo ur selection criteria (e.g. Wig- 
gleZ; iGlazebrook et al] l2007t ) or require a bright emission 
line to obtain a redshift (e.g. the Joint Dark Energy Mission, 
JDEM3) this is not an option. In this section, we will first 
relate the distribution of disk orientations to b K and then 
derive the A — b K relation for general inclination-dependent 
flux 4>(i). We will then consider a few toy models of 
and their implications. More complex models have been con- 
sidered else where (e.g. iDisnev et al.|[l989l ; iGiovanelli et al] 
Il994l,ll995h . 

We consider a distribution of disks whose normal vec- 
tors L have probability 

P(t) = -^(l + BB ij iitt), (55) 

where B describes the extent to which large-scale tidal fields 
affect the orientation of galactic disks. The anisotropy of the 
form Eq. ([55} is the only correction that is allowed by sym- 
metry to first order in the large-scale tidal field. Such an 

7 URL: http://jdem.gsfc.nasa.gov/ 
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anisotropic alignment can occur due to no nlinear evolution 
in tidal torque models l|Hui fc Zhandl2008l ). For a geometri- 
cally thin axisymmetric disk perpendicular to L, the incli- 
nation is cosi = I/3 and the axis ratio is b/a = cosi. Thus 
the ellipticity is 



iP/af 



1 + (&/a) 2 



L\ + L\ 
1 + Ll ' 



(56) 



The position angle <j> of the apparent major axis is 7r/2 away 
from the position angle of the projected angular momentum 
vector, tan _1 (I/2/I/i). Thus the components e\ — ecos2<^> 
and e 2 = e sin 2(f) of the ellipticity tensor are 



ei 



and e 2 



2L X L 2 



(57) 



The mean ellipticity can be obtained by integrating the 
probability distribution Eq. (|55|l for L over the unit sphere. 
That is, 



(ei + ie 2 ) = - 



1 



f s2 p(i)d*L 



f ik 

Js? 1 



+ iL 2 



-(l + BsyLiL f )p(i)d 2 L, 



(58) 

where p(i) is the selection probability for a galaxy at incli- 
nation i. We switch to spherical polar coordinates (i, (j>), and 
performing the cf> integration we obtain 



(ei) = -K(su - S22); (e 2 ) = ~2Ksi2, 



where 



K 



Jq sin 5 i (1 + cos 2 i) 1 p(i) di 
4 fg p(i) sin i di 



(59) 



(60) 



If the anisotropic selection is a small effect so that 
p(i) ^constant, then we have K = tt/4 - 2/3 « 0.1190 
In models with large anisotropic selection effects, K might 
be decreased - for example, in the extreme example below of 
an optically thick uniform disk whose apparent luminosity is 
proportional to cos i, if the slope of the luminosity function is 
r\ = d In n/d In F m i n = 2 then p(i) oc cos 2 i. This would result 
in K = 0.044, nearly a factor of 3 smaller. We will param- 
eterize this uncertainty with the parameter = 7f/0.119, 
which is equal to unity for isotropically selected thin disks 
and declines if face-on disks are preferentially selected[f| Note 
also that the effective K would be decreased if we took into 
account that even a perfectly edge-on disk does not have 
e = 1 because of its finite thickness and the contribution 
of the bulge. Both effects imply that in practice one should 
have K a < 1. 

For a Fourier mode k |] ei in the plane of the sky, we 
find ei(fe) = KBSm(k), and hence 



KB 



-0.068.KaB. 



(61) 



cos i and decomposing in 



8 This is obtained by substituting x 
partial fractions. 

9 By inserting <5- functions in p(i), we see that — > if only 
face-on galaxies are selected, and K a -> (n - 8/3) -1 2.11 if 
only edge-on disks are selected. 



This establishes the rel ation between th e disk orientations 
and the b K parameter of Be rnsteinl l|2008l ). 

In order to proceed, we must determine how the selec- 
tion probability e is related to B. We assume a flux function 
with slope dlnn/dlnFinin = —r\. Then the number density 
of galaxies per logarithmic range in intrinsic flux F\ per unit 
solid angle of orientation is 



N(Fi,L) oc Ff (1 + BsijLiL 



(62) 



The observed flux is F = fi$(i). The number density of 
galaxies above some threshold flux Fq is then 



N(> Fo) oc / d 2 L 



dlnFiFf" [1 + BsijLiL. 



(63) 

One can separate the integral over d 2 i into an integral over 
inclination sini = L3 and position angle <j>- The position 
angle integrates out trivially, 

N(> Fo) oc / sinidi / dlnFFr 11 

JO JF a /i>(i) 



1 + Bs 33 cos 2 i + -^(sn + S22) sin 2 i 



(64) 



Using the tracelessness of Sij to simplify the last term, and 
performing the Fi integral, we get 

N(>F )(x mi)]" [I+BS33P2 (cos i)] sini di. (65) 



The anisotropic dependence is that due to S33: 



r 7r [$(i)] I) P 2 (cosi)sinidi 
£ = Bs 33- H — TWT^T^Z. ■ ■ ■ ( 66 ) 



Io[^( i )] v sinidi 



Thus, defining 



Jo[$(i)] r >P2(cosi)sinidi 

^ = — — r^nrrPvu. — > ( 67 ) 



So [®{i)] v sinidi 
we see from Eq. ((28)) that A = —tpB or 



A 



(68) 



Since P2(cosi) is bounded between — | and 1, it follows that 
— 5 ^ V 1 ^ 1; f° r selection in favour of face-on galaxies 



Tp > 0. The parameter is obtained from Eq. (|60[) using 
the selection probability, which is p(i) oc [<I>(i)]''. 

In order to use Eq. (|68[) we need to construct a model for 
the angular distribution of emitted radiation <3?(i). We con- 
sider several toy examples. These are simple plane-parallel 
models in which we consider the dust optical depth r 
through the galactic disk. This need not equal the optical 
depth relevant for extinction-correcting the stellar contin- 
uum or emission lines. For example, in an [O 11] or Ha flux- 
limited survey where there may be significant internal ex- 
tinction in H n regions, the analysis below should include 
only the diffuse contribution due to propagation from the 
disk (which is affected by disk inclination), rather than the 
"total" optical depth. 

Uniform slab: We suppose the emitters have the same 
vertical distribution as the diffuse dust, and the vertical op- 
tical depth through the disk is r. In this case the observed 
optical depth at inclination i is r|seci|. Then &(i) is pro- 
portional to the probability of a photon emitted at random 
location in the disk reaching us without being absorbed, i.e. 
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$(«) = (l-e" T|scc!| )/(r|seci|). The integrals in Eq. (JBTJl are 
not represented by any commonly used function and must 
be evaluated numerically, except in the limiting case of an 
optically thick slab (r ^> 1). 

Uniform, optically thick slab: If the disk is optically 
thick to dust absorption and has emitters (young stars or H 
II regions) with the same vertical distribution as the diffuse 
dust, then the surface brightness of the disk is inclination- 
independent. This gives a flux proportional to projected area 
or <f>(i) oc | cos i| . Both the numerator and denominator in 
Eq. (|67|l are then polynomials in cos i; the integral evaluates 
in closed form to 



n 



(69) 



r/ + 3 

There is no closed-form analytic expression for K a for gen- 
eral r\. However for integer r\ one may use the substitution 
x = cos i and integrate via partial fractions to get 



K=l 



77 

4 


2 
3 


'/ 


= 


In 2 


5 


V 


= i 


12 

5 

17 
12 




V 


= 2 


2 In 2 


V 


= 3 


!*- 


82 
21 


n 


= 4 



(70) 



which imply K & = 1, 0.574, 0.369, 0.256, and 0.187, respec- 
tively. 

Emitting sheet embedded in absorbing slab: If the disk 
has a thin emitting layer sandwiched between two ab- 
sorbing slabs each of optical depth r/2, then <E>(i) oc 
With the substitution of x = cosi, Eq. (|67[) 



T SCC I 



1/2 



is transformed into an exp onential integral (Eq. 5.1.4 of 
lAbramowitz fc Stegun| [l972). giving 



3E 4 ( V t/2) 
2E 2 (r,r/2) 



-\ryr ln(2.42r?r) 



T)T < 1 
■qr > 1 



(71) 

The asymptotic expans i ons u se Eqs. (5.1.12) and (5.1.52) of 

lAbramowitz fc Stegunl (ll972T). 

Power-law model: iGiovanelli et all (|l994T ) estimated 
that the total (exponential model) magnitude of local Sc 
galaxies suffers an inclination-dependent correction AM = 
— 7log 10 (a/6) in the I band, where 7 = 1.02±0.08 was found 
by minimizing the x 2 of the corrected Tully-Fisher relation 
fit. This implies <E>(z) oc | cos' 

i> = 



o.4 7 , 



and 



0.4?77 



O.4777 + 3 



(72) 



Note that in all of these cases additional corrections 
apply if isophotal or aper ture magnitudes are u sed in place 
of exponential model fits (|Giovanelli et al .11 19941 ). 

We consider a slope of the cumulative luminosity func- 
tion 77 = — d In n/d In f m i n « 2. This is roughly appropri- 
ate for a JDEM-type survey using the Ha em ission line 
to a density of several x 10 ~ 4 Mpc -3 at z ~ 1 jYan et al.l 
19991; ISumiovshi et alj|2009ft and z ~ 2 (|Reddv et al.ll2008l ; 



Geach et al.l 20081 ). In this case, the optically thick slab 
model would predict ip = 0.40 and K a = 0.37, or if)/K a = 
1.1. The sheet-in-slab model predicts ^ = 0.10 and 7f a = 0.8 
or ifc/K* = 0.12 at r = 0.1, rising to i>/K a = 0.5 at r = 0.5. 
More values are plotted in Figure[3] For the power-law model 
with 7 £3 1.0 appropria te for local Sc galaxies in the I band 
l|Giovanelli et al.ll 19941) . we find ip = 0.21 and A' a = 0.63, or 
ip/Ka = 0.33; a somewhat larger value might be expected 



Anisotropic selection factors for disk galaxy models 




Figure 3. The orientation-dependent selection factor ip/Ka, for 
both the uniform slab and sheet-in-slab models for r\ = 2, as a 
function of disk optical depth. 



at bluer wavelengths since there should be more internal ex- 
tinction in the galaxy, and one might expect younger stars 
and Hll regions to have a smaller scale height. It is readily 
apparent that many models with significant optical depth 
give ip/Ka in the range of a few tenths. 

The last ingredient needed to make a prediction for A is 
a prediction for or measurement of b K . Unlike the case of el- 
liptical galaxies, for disk galaxies there is at present no clear 
consensus on whether b K 7^ 0. Theoretically, b K is exactly 
zero to lowest order in tidal torque theory because a collaps- 
ing protogalaxy requires not just a tidal field to spin up but 
also an anisotropic quadrupole moment of the mass distri- 
bution; the combination leads to alignments that do not cor- 
relate with the densit y perturbations and go rapidly to zero 
in th e linear regime ([Crittenden et al.l l200ll; ICatelan et al.l 
200 it iMackey et al.ll2002l: Iffirata fc Seliakll2004l ). However. 



Hui fc Zhang (|2008) showed that when nonlinear evolution 



is considered, there is a squeezed bispectrum configuration 
of small-scale tidal fields, quadrupole moments of collapsing 
structures, and linear-scale density perturbations, so that 
tidal torque theory predicts b K > 0. iHevmans et alj (|2006l ) 
used A^-body simulations to search for correlations of spin 
axes with the density (technically the cosmic shear) field and 
found none. 

On the observational side, some early works found evi- 
dence for disk galaxy alignments with large-scale structure 
i Lee fc Penll2002l ; lNavarro et al.ll2004) and interpreted their 
results as supportive of tidal torque theory. On the other 
hand, studies with the much larger SDSS samples of blue 
galaxies have turned up no evi d ence for density-ell i pticit y 
correlations l|Hirata et al.l 120071 ; iFaltenbacher et al.M2008l ). 
even w hile confirmin g previous detections of LRG align- 
ments llBinggelil ll982h at high significance. There have also 
been several claimed detections of disk galaxy orientations 
with the intermediate axis of the tidal field o r with void 
shells, as would occur in tidal torque theory l|Lee fc Pen! 
l2002tlTruiiUo et al.ll2006l ; iLee fc Erdogdull20071 ) . These were 
all based on 3-dimensional modeling of the tidal field in 
which the line of sight is a preferred axis for both the tidal 
field (due to redshift space effects) and the galaxy angu- 
lar momentum reconstruction; for this reason it is desirable 
to confirm these results with statistics such as position an- 
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gle co rrelations where n o direction is preferred by the algo- 
rithm. iLee fc Pen! (|2007l ) searched for position angle correla- 
tions of neighbouring blue galaxies in SDSS; they observed 
a 3cr positive correlation in their innermost radial bin (3D 
redshift-space sepa ration s < 3ft~ M pc) with nondetection 
at s > 3ft" 1 Mpc. ILee fc Pen! (|2007t) noted that this scale 
dependence is consistent with a quadratic tidal alignment 
model, although the statistical significance is marginal. This 
would lead to a next-order term involving SijSki in Eq. (I28p . 
whose correlations drop off rapidly at large scales. This ad- 
ditional term would not affect RSDs in the large-scale limit, 
but the quasilinear regime would be affected. 

Faced with these uncertainties, we consider as a worst- 
case scenario for the linear b K coefficient for disk galaxy 
alignm ents the "pessimistic" models from iHirata et al.l 
((2003). The pessimistic model corresponds to the 2o up- 
per limit on alignments of ~ L* late-type galaxies; this 
is \w s+ {r p = 20ft" 1 Mpc, z = 0.3) | = 0.028ft" 1 Mpc, or 
\b K \ = 0.008. Assuming ip/K a ~ 0.33 in accordance with 
the power-law (7 = f) model, this leads to \A\ ~ 0.039 and 
a fractional contamination to the redshift-space distortion 
|A|// of 6 per cent. This should not be taken to mean that 
there is contamination at this level, but rather as an ex- 
ample that may be appropriate if disk galaxy alignments 
are near present upper limits. The true value of |^4|// may 
be much smaller (although larger values at high redshift, or 
larger ip/K a for some types of selection, are possible). 



6 REDSHIFT-SPACE DISTORTIONS VERSUS 
COSMIC SHEAR 

ft is of interest to compare the effect of tidal galaxy align- 
ments on RSDs to their better-known effect on CS measure- 
ments, assuming that one uses similar samples of galaxies. 
We will conclude that for the disk galaxies, even if ip/K a is of 
order unity, then for low-redshift measurements the cosmic 
shear is more affected than the RSDs. However at redshifts 
of order unity, the fractional contamination to the RSDs is 
generally ip/Ka, times the fractional contamination to the 
cosmic shear. Since we have seen that for many simple disk 
models ip/K^ is a few tenths, and since future galaxy surveys 
designed to probe dark energy will be exploring the z ~ f 
range, tidal alignments should if uncorrected be considered 
a similar source of error for both types of measurements. 

The contamination of the RSD measurement is very 
simply 

A(z) ^ A(z) 



A\n[f(z)G(z)] = 



f( Z ) [On(z)P 



(73) 



The contamination of the CS signal is more compli- 
cated. We consider the shear power spectrum for a broad 
redshift distribution, which for the pure cosmic shear case 
is given by 

MOT; 



Up — 



I 



k= - ; To ■ 
r 



dr. 



(74) 



where P m is the matter power spectrum at the wavenumber 
k and conformal time to — r, the conformal time today is to, 
the lensing window function is 



W(r) = -Sl m H%[l + z(r)]r 



-n(r s )dr s , (75) 



n(r s ) is the probability distribution for sources to lie at 
distance r a and we have assumed a flat universe (e.g. 
IHirata fc SeliaJd 12004 ) . Now suppose the galax ie s hav e an 
intrinsic alignment described by the iBernsteinl (|2008l ) pa- 
rameter b K and (as appropriate for large scales where the in- 
trinsic alignments are due to the local tidal field) P s ^i (fc) = 
b K P m (k). Then the leading-o rder contamination in b K is the 
"GI" or interference term of IHirata fc Seliakl l|2004 l: 



cr = 2 



^p^ bK pJk=L,T0 

r 4 \ r 



dr. 



(76) 



The evaluation of the fractional contamination to the 
lensing power spectrum C GI /C GG requires a numerical com- 
putation in general. However for conceptual purposes it 
helps to take the limit of nearby sources, i.e. z < 1 so that 
r ss H ~ 1 z, and take P m (ft) to be a power law P m (fc) = 

off 

Nk n * . For definiteness we take a source distribution 



r» — r / r 



n(rs) = 27f 

which implies a lensing window function 

WM - -<) 11-, ■ I — - > \ ~ r / r * 



(77) 



(78) 



The pure lensing contribution to the power spectrum is 
16 



C 



: Q 2 m H^rl- nt 



u (u + 2) e 



dtt 
(79) 



with the substitution u = r/r*; this simplifies to 

T 



C, 



GG 



—nLHoMn 

16 

r(3 -nf) 



r 



off 2 15 off ... „ 

-n s - —n s + 13 



23-nJ ff 

Finally the GI contribution is 



Go — —£l m H, 



(80) 



b K Nr* & I 3 / u (u + 2)e u di 



which simplifies to 
3 



G 



^HlbMr: 1 ^ t 



. r(4 - < 



1-inf 



The ratio is 

cf 



4(3 -nf )(8-nf ) b K 

Cf G 3(52 - 15n| ff + nf 2 ) f2 m # 2 r2 



(81) 



(82) 



(83) 



The prefactor is a combination of dimensionless integrals 
that depends on n^ s ; for typical values of the spectral index 
in the quasilinear regime nf ~ —1.3, it is 0.73. Also the 
median redshift is z mc d ~ 2.67// ~ 1 r*. This implies 
-<gi u 



C 



5.2 



Q 7 2 ' 



(84) 



The fractional error in as depends on how Ce scales with ag. 
In the linear regime, C GG oc cr|, but in practical cases, CS 
uses the quasilinear regime and the actual scaling is closer 
to Cf G oc <rf. Thus if one uses cosmic shear to normalize 
G(z -C 1), then the inferred error on the normalization of 
the power spectrum is ~ C GI /3C GG , 

Aln G^1.7 ?r ^. (85) 
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The ratio of the contamination in the RSD measurement to 
that in the CS measurement is then, at low redshifts, 



Aln(/G)| R sp 
AlnGlcs 



n « o - 4 2 
U.O — Si m Z mcd . 

b K 



(86) 



(The factor of 0.6 depends on the shape of the redshift dis- 
tribution and should be treated only as a rough estimate.) 

In the disk galaxy case, we then have from Eq. (|68[1 that 
0.07A/b K « V/A'a, so 



Aln(/G)l RS i 
AlnGlcs 



o j> pO.4.2 

^ m 2 mccl • 



(87) 



Thus at low redshift z < 1 and ip/K a ^ 1 the RSD measure- 
ment suffers less contamination than the CS measurement. 
This is because the lensing strength for low-redshift surveys 
is small, and hence the ratio of intrinsic galaxy ellipticities 
to lensing is enhanced. However, when the source redshift 
becomes of order unity, Eq. (|86[) implies that if ip/K a is of 
order unity, then the RSD and CS measurements suffer sim- 
ilar levels of contamination. This statement should be inter- 
preted only in an order-of-magnitude sense, since Eq. (|8T[) 
was derived by leaving out factors of 1 + z; more detailed 
comparisons would require knowledge of how b K , xp/K^, etc. 
vary with redshift. Nevertheless, we can conclude that if we 
use disk galaxies at z ~ 0(1), and if the anisotropic selec- 
tion parameter ip/K a ~ 0(1), the fractional contamination 
of the redshift- space distortion measurement due to galaxy 
alignments is inherently of the same order of magnitude as 
the fractional contamination of the cosmic shear measure- 
ment. 

The simple analysis presented here only refers to the un- 
corrected magnitude of the systematic introduced by tidal 
alignments. It is much harder to compare the impact of tidal 
alignments on RSD and CS measurements after corrections 
are applied, because the nature of the correction is different 
in the two cases. In CS, intrinsic alignments, even in the 
nonlinear regime, can be ge ometrically projected out based 
on the redshift dependence (|Hirata fc Se liak 20041), albeit at 
the expense of increased error bars an d tighter requirement s 
on photometric redshift performance (|Bridle fc Kindl2007h . 
In contrast, systematics control for RSD would have to in- 
volve a combination of choosing galaxies that are weakly 
aligned (6 K ~ 0), have selection criteria that are close to 
isotropic (x, ip ~ 0), and estimating the residual contamina- 
tion as well as possible. It would be premature at present to 
conclude that one of these strategies will be more successful 
than the other. Moreover, since CS and RSD do not measure 
the same function (CS measures G whereas RSD measures 
/G) they are not necessarily competitors. 



surveys using the brightest LRGs for some types of selec- 
tion criteria the effect could be ~ 10 per cent, whereas 
for nearby disk galaxies using the 2a upper bounds on 
tidal alignments and making reasonable assumptions about 
inclination-dependent selection effects we expect ~ 6 per 
cent. Realistic models would probably give smaller contami- 
nation, since the tidal alignment effect decreases as one goes 
down the LRG luminosity function and our result for disk 
galaxies uses an upper limit. The effect is probably unim- 
portant for current RSD measurements, but could represent 
a problem for future large-scale structure surveys. 

There are several strategies to reduce contamination 
of the RSD by tidal alignments. One could choose sam- 
ples of galaxies with weaker alignments; in LRG surveys 
this means going to fainter galaxies, whereas for late-type 
galaxies we do not have a detection of tidal alignments so no 
specific subsample can be identified as "best" at this time. 
Secondly, one could try to reduce anisotropic selection ef- 
fects; for optically thin galaxies such as (most) LRGs, this 
would mean using model magnitudes for selection. For late- 
type galaxies whose optical luminosities are significantly af- 
fecte d by dust this is harder, bu t to some extent is possible 
(e.g. lUnterborn fc Rvdenll20o"^ ). A final strategy would be 
to take observed density-ellipticity correlation functions and 
attempt a model-dependent calculation of the A parameter 
by methods similar to those used in this paper (but more 
sophisticated) . 

In this paper we have focused our attention on the ef- 
fect of tidal alignments in the linear regime. The same pro- 
cesses occur in the quasilinear and nonlinear regimes, al- 
though there it is even more difficult to make predictions. 
We anticipate that future RSD work will make extensive 
use of the quasilinear regime as a test for systematics (e.g. 
improper finger-of-God compression, second-order biasing, 
etc.), and possibly to correct "linear" scales for such effects. 
The implications of tidal alignments for quasilinear RSDs 
are left to future work. 
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7 DISCUSSION 

In this paper, we showed that tidal alignments of galaxies 
can be a source of systematic error for the redshift-space 
distortion technique to measure the growth of cosmic struc- 
tures. This is because a combination of tidal alignments and 
a viewing-direction-dependent selection function results in 
a different observed power spectrum for Fourier modes in 
the radial and transverse directions. We have considered 
several toy models for the effect. The amount of contam- 
ination is highly uncertain and survey-dependent; but for 
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APPENDIX A: REALITY CONDITIONS 

This appendix derives the reality condition satisfied by the 
power spectra in Eq. (|18|l . We denote by R the rotation 
matrix corresponding to — fc with — Rfc = ez introduced in 
Section f3. 21 and define S = RR _1 . 

The coefficients in Eq. JT7} are related since e(h\x) is 
real. This implies e(njfc) = e*(n| — fc) and so 



E 



(fc)y im (Rn) 

v / 2TTT 



E 



i'e7 m (-fc)yr m (Rn) 

v / 2TTT 



Im v I m 

Defining u) = Rn, we convert Eq. (|A1|) into 



E 



i l ei m (k)Y lm (u) 



1 , 



gjm(-fc)yjm(S^) 

V / 2TTT 



(Al) 



(A2) 



By definition —Rfc = e$ and Rfc = £3, so Se3 = — £3. 
Since S GSO(3) and has a —1 eigenvalue, S must be a rota- 
tion by 7r radians around an axis orthogonal to £3. We sup- 
pose this axis is at position angle <E>, i.e. is cos$ei +sin«I>e2. 
In spherical coordinates, S therefore takes the point uj = 
(6, <f)) and maps it to the point = (n — 9, 2$ — <f>). We 
conclude that 

YT m (Sw) = YZ m (4>-6,2&-4>) 
= Y lim (4>-6,-2& + 4>) 
= (-l) i Y ! , m (0 j7 r-2$ + <£) 

= (-l) i+m e- 2im *Y ijm ( W ). (A3) 

Here the second line used the Condon-Shortley phase con- 
vention, the third used the ( — 1)' parity of the spherical har- 
monics, and the fourth used the (^-dependence. Comparing 
to Eq. (|A2[) . we obtain 



ezm(fc) = (-1) 



m — 2im"i> 



(A4) 



Thus the fc and — fc components are related, up to a phase 
$ that depends on the choice of R and R. We may then take 
the expectation values of products of Fourier modes, 

/—* /I \ - /I \\ / -i\m — m 2i(m — m )(& 

(e im (fc)e,/ m /(fc)} = (-1) e k ' 

x(e* m (-fe)ei'm'(- fe ')>*- ( A5 ) 

Since by rotational invariance the power spectra are only 
nonzero for m = m', and the power spectra must be the 
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same for k and —k, comparison to Eq. (|18[1 shows that the 
power spectra P £ m (k) must be real. 



APPENDIX B: GENERAL INTRINSIC 
ALIGNMENT POWER SPECTRUM MODELS 



This appendix generalizes the results of Sec. !3.3l to intrinsic 
alignment models with non-scalar contributions (m 7^ as 
well as m — 0) and galaxy models lacking inversion symme- 
try (odd as well as even I). 

Equation (|23|l generalizes to 



P S (k) 



(6 + /p 2 ) 2 P m (fc) 

OO 

+2(& + / M 2 )5R]T(- 

1 

.1-1' 47rPj 



1=1 

-yll m 



Pi e (fc)y !0 (Re 3 ) 



21 + 1 
(fc)y,: n (Re 3 )^ m (Re 3 



X /(2Z + 1)(2«' + 1) 



(Bl) 



The cross-term in Eq. (IB1|I can be simplified by replacing 
Yio with a Legendre polynomial using Eq. (J25J . The real 
part K kills the odd-i terms. Physically the reason these 
do not contribute to the observed galaxy power spectrum 
(even if present) is that the odd-Z intrinsic alignments are 
7r/2 radians out of phase with the matter perturbations, 
hence there is no interference in the power spectrum. 

The pure intrinsic alignment term in Eq. (|B1|) can also 
be simplified, albeit with more work. The I — I' odd terms 
cancel because (1,1') and enter with opposite sign 

(again a manifestation of them being n/2 radians out of 
phase) . The I — I' even terms simplify using 



47ryi^(Re 3 )yrm(Re 3 
^/(2/ + l)(2Z' + l) 



= Nt m Nt, m pr(j*)pr(j*), (B2) 



where Ni m = y(T-j- m)\/(l — m)\. Thus we finally have 
P g (fc) = (& + /p 2 ) 2 P m (fc) 

+2(& + / M 2 ) y. (-v l/2p Uk)Pi(») 



+ 



E 



l ^ 2 , even 

(_l)y-^)/2p«'0 (A;) p i(M) p i;(At) 



IV ,l — V even 



+2E E (-i) (( - ! ' J 

m>0 W ,1—1' even 



N lm N Vm Pl l ' m {k) 



(B3) 



Here we have split off the m = (scalar) pure intrinsic 
alignment term, and in the m 7^ terms we have used the 
equality of ±m contributions to restrict the sum to m > 
and multiply by 2. 



We begin by defining a new matrix 

u = |Wi2 ivy 22 v^ 32 
V iw 33 



(Cl) 



Then U + U T = W and to first order in W, 

p 2 = s-(l + W) _1 s w s-(l-U-U T )s^ s-(l-U) T (l-U)s. 

(C2) 

so 



|(I-U) S | 



(C3) 



We now define the vector c = (I — U)sx, which lies in the 
12-plane. Equation (|36p simplifies to 



J(\c + a 3 (\-U)e 3 \) ds 3 , 



(C4) 



This would be the projected image of the spherical galaxy 
model (i.e. with W = 0) at position c if it were not for 
the presence of U. To first order, the U13 and I7 23 terms 
produce contributions to the integral odd in s 3 so they can 
be dropped, leading to 



I(s±) = J J (|c + S3 (l - (7 33 )e 3 |) ds 3 



(C5) 



The factor of 1 — C/ 33 can be eliminate by rescaling the radial 
variable to (1 — i7 33 )s 3 , thereby introducing a factor of the 
Jacobian (1 - t/ 33 ) _1 = 1 + \ W^. This leads to Eq. J37|. 



APPENDIX D: 
GALAXY 



ALIGNMENT OF TRIAXIAL 



In Sec. 15. II we described the triaxiality of an elliptical galaxy 
by a traceless-symmetric matrix W. It was argued that if 
there exists a valid Taylor expansion of W in terms of the 
local tidal field Sy, then one may write: 



(Wij) — 2Bsij + higer order terms, 



(Dl) 



i.e. Eq. (|50p . The purpose of this appendix is to provide 
an alternative derivation based on the distribution of Euler 
angles of the galaxy. 

Since W is a traceless-symmetric matrix, it can be de- 
scribed by its diagonalized form, 



W 



Q T AQ, 



(D2) 



where 



APPENDIX C: PROJECTED ELLIPTICAL 
GALAXY 

This appendix derives the projected elliptical galaxy image 
I(s±), Eq. (1371) . The solution presented here is equivalent 
to that of IStarkl (|l977r ) in the limit of small W (t, U — » 1 in 
IStarkHl977T notation). However for clarity we present it here 
in a different coordinate system and using notation that is 
well-suited for the tidal alignment problem. 



A = 



Ai 











A 2 











A3 



(D3) 



is diagonal and Q 6 SO (3) is an orthogonal matrix specifying 
the orientation of the galaxy. Specifically, Qin rotates from 
the "lab" (unprimed) frame to a coordinate system aligned 
with the three principal axes of the galaxy (primed). It can 
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be expressed in terms of its Euler angles: 

Qii = — sin ip sin <f> + cos 6 cos ip cos <p, 

Q12 = sin')/) cos cp + cos 6 cos ip sin <p, 

Q13 = — sin # cos V", 

Q21 = — cos ip sin — cos 6 sin ?/> cos 0, 

Q22 = cos V> cos <f) — cos 6> sin ip sin 0, 

Q23 = sinSsini/), 

Q31 = sin # cos <f>, 

Q32 = sin 6 sin 0, and 

Q33 = cos#. 



(D4) 



Tracelessness guarantees Ai + A2 + A3 = and we can al- 
ways choose Ai A2 ^ A3 An oblate galaxy would have 
Ai < A2 = A3 and a prolate galaxy would have Ai = A2 < A3. 
We consider the case of a population of galaxies with fixed 
3D axis ratios, i.e. fixed {A;/}. The generalization to a dis- 
tribution of axis ratios is immediate. 

In the absence of a tidal field, the orientation is 
isotropic, i.e. the probability distribution for the Euler an- 
gles is simply the group volume element divided by the vol- 



dP = 



> A6 dtp dcj) 
"8^2 ' 



(D5) 



where as usual the volume is taken over the range ^ 6 ^ n, 
^ ip < 2-7T, and ^ <p < 2n. In the presence of a tidal field, 
a general linear-order correction in Sij is: 



dP= [i + sijdjie,^)] 



sin 6 d9 dtp d<p 
8^ ' 



(D6) 



The specific form of the function Cy can be constrained by 
rotational invariance arguments. The product SijCij{9, <p, ip) 
can depend only on the components of Sij in the galaxy 
frame, which are s^/ = QiuQy jSij . That is, 



sin 6 d<9 dip dcp 
8^ 



dP = [1 + Q iH {6, <j>, iP)Q fj (6, <j>, 

(D7) 

where now all of the dependence on the Euler angles is 
packaged into Q, and CVy has no dependence on the Eu- 
ler parameters. The matrix CVy is contracted with the 
traceless-symmetric tensor Sj/y and hence can be taken to 
be traceless-symmetric without loss of generality. For galax- 
ies with triaxial symmetry, the 180° rotation symmetries 
around each axis guarantee that CVy will be diagonal 1"1 
Hence C is traceless-diagonal and has 2 degrees of freedom. 
Thus 2 parameters are required to fully describe the linear 
alignment of a triaxial object by a tidal field. 

(In the special case of an axisymmetric galaxy, e.g. a 
prolate galaxy with Ai = A2, we have Cu = C22 = — ^33. 
In this case, there is only 1 parameter required to describe 
C and hence specify the intrinsic alignment model.) 

We next return to the problem of finding the average 



10 Since the eigenvectors are only defined up to a sign, there are 
4 equivalent choices of Q (not 2 3 = 8 since we restrict to proper 
rotations). In what follows we assign one of these 4 choices to 
each galaxy with equal probability. 

11 e.g. C12 flips sign under 180° rotation around the galaxy 1- 
axis. 



value of W: 

(Wij) = Ai'j' (Qi'iQj'j) 



h-i'jl I Qi'iQj'j[^- + Qk'kQl'l 3 klCk'l' 



x - 



sin 9 d9 dtp dcp 



8tt 2 



iji'j' 1 J 

where we have defined 



(2) 



n 



and similarly for II 



Qi'iQ 



(4) 

ijkli' j' k' ' V ' 



ijkli 1 j'k'l'-^-i' 'f ' S klCk>V , 

sin 9 d6 dtp dtp 



8tt 2 



We next need to evaluate n^, j( and ^-fjlwi'k'l' 

t(2) 



(D8) 
(D9) 
Con- 



sider first .;. Since the group volume element is invari- 
ant under right multiplication, it follows that for any rota- 
tion matrix O we may replace Q — > QO in Eq. ()D9|I . thus 
establishing: 



n (2) , ., = O ia O lb Il 



(2) 

abi' j' 



(D10) 



Now we know that any 3x3 matrix X that is rotationally 
invariant, i.e. X = OXO T for all O GSO(3), is proportional 



to the identity. For fixed i'j', H^},., forms such a 3 x 3 ma- 

trix and hence 11^,., = Pi/jiSij for some Pi'ji. A similar 
argument using left multiplication shows that p^ji is pro- 
portional to the identity, so 



II 



(2) 



cSi/jiSi 



(Dll) 



Direct integration of Eq. (|D4f) for II3I33 = | enables us to 



solve for c = | : 



II 



(2) 



-Si'j'Si 



The similar argument for 'n^'j e ii'j'k'i' is more compli- 
cated. There are three linearly independent fourth-rank ten- 
sor Xijki that are invariant under rotations, i.e. X^m = 
OiaOjbOkcOidXatcd- A general such tensor can be written 



(D12) 



Xijki = aiSijSki + a>2Sik5ji + a-jSuSjk- 
Thus the analogous equation to Eq. (|D 11 f) is 



(D13) 



n 



(4) 

ijkli' j' k' I' 



Ci8ij5 k i8i'j>5k'i' + C25ik8ji5i'ji8 k ii, 
+C38u8jk5i'j'5k'i' + C4,8ij8 k i8i'k'8j'i' 
+Cs8ik8ji8 i iki8yii + C(,8u8jk8i' k '8jii' 
+c 7 8 i j8kiS i ii l 8 j i k i + cg8i k 8ji8i>i'8ji k > 
+cg8a8j k 8iii'8ji k i . (D14) 

This can be simplified using the permutation symmetries 
such as II.*: , ., .,, ,„ 

zj Kit J k t 



IlS,.,,,^,,. These imply ci = c 5 = c 9 



and C2 



n 



(4) 



C'3 



C4 



jiklj' 
C6 = C7 



C8. With the constraints 



33333333 = I = 3ci + 6c 2 and n^Mnsa = = ci, we then 
find ci = and C2 



30 ' 



12 This foll ows from the general theorem, e.g. Theorem 2. 9. A of 
IWevll lll939l '). that all rotation-invariant tensors can be expressed 
as polynomials in terms of the metric tensor Sij and the determi- 
nant tensor 6y&. 
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16 Hirata 



Combined with the tracelessness of A and s, diagonality 
of A, and symmetry of s, this implies 



The quantity in parentheses is then merely a constant, thus 
verifying Eq. (f50)l . 

It is interesting to note that to lowest order in the devi- 
ation from sphericity, the viewing direction-dependent selec- 
tion function will depend only on (Wij). Thus the contami- 
nation to RSD (and to CS, since to lowest order the intrinsic 
ellipticity of a galaxy is Wn — W22 + 2iWi,2) depends only 
on 



Therefore, even though the intrinsic alignments of a triaxial 
object are described by 2 parameters (the linearly indepen- 
dent values of CW), only 1 linear combination of these is 
relevant to our analysis. 




(D15) 




(D16) 
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